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Abstract 

We find a two-degree-of-freedom Hamiltonian for the time-symmetric problem of straight line 
motion of two electrons in direct relativistic interaction. This time-symmetric dynamical system 
appeared 100 years ago and it was popularized in the 1940's by the work of Wheeler and Feynman 
in electrodynamics, which was left incomplete due to the lack of a Hamiltonian description. The 
form of our Hamiltonian is such that the action of a Lorentz transformation is explicitly described 
by a canonical transformation (with rescaling of the evolution parameter). The method is closed 
and defines the Hamitonian in implicit form without power expansions. We outline the method 
with an emphasis on the physics of this complex conservative dynamical system. The Hamiltonian 
orbits are calculated numerically at low energies using a self-consistent steepest-descent method (a 
stable numerical method that chooses only the nonrunaway solution). The two-degree-of-freedom 
Hamiltonian suggests a simple prescription for the canonical quantization of the relativistic two- 
body problem. 
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I. INTRODUCTION 



The class of equivariant dynamical systems under the Poincare group has enormous rele- 
vance to physics and yet, to date, only the one-body relativistic motion is fully understood. 
Already with two bodies in relativistic motion, one encounters the no-interaction theorem: a 

Boup theoretical obstacle to the Hamiltonian description of relativistic two-particle motion 
. The no- interaction theorem can be overcome by covariant constraint dynamics |2|, but 
one is left with the few cases where the constraint scheme closes. For example, for an equiv- 
ariant physical theory like the time-symmetric electrodynamics [3], a constraint description 
is unknown. In this paper we present a reduction of the time-symmetric two-body problem 
of the relativistic action-at-a-distance electrodynamics to a two-degree-of-freedom Hamilto- 
nian system along the nonrunaway solutions. The form of the Hamiltonian is such that a 
Lorentz transformation is explicitly described by a canonical transformation with rescaling 
of the evolution parameter. The Hamiltonian orbits are calculated numerically by a numer- 
ically stable self-consistent method that uses steepest-descent quenching and chooses only 
the nonrunaway orbits. 

In 1903, Schwarzchild proposed a relativistic type of interaction between charges that 
was time reversible precisely because it involved retarded and advanced interactions sym- 
metrically ^J. The same model reappeared in the 1920s in the work of Tetrode and Fokker 
and it finally became an interesting physical theory after Wheeler and Feynman showed 
that this direct-interaction theory can describe all the classical electromagnetic phenomena 
(i.e. the classical laws of Coulomb, Faraday, Ampere, and Biot-Savart) P,!^- Another ac- 
complishment was that Wheeler and Feynman showed in 1945 that in a certain limit where 
the electron practically interacts with a completely absorbing universe, the response of this 
universe to the electron's field is equivalent to the local Lorentz-Dirac self-interaction theory 
0] without the need of mass renormalization 3] . It is amusing to understand that the clas- 
sical radiative phenomena of Maxwell's electrodynamics can be described as a limiting case 
of this direct-interaction theory (complete absorption is added to the theory as a simple 
model to uncouple it from the detailed neutral-delay dynamics of the other charges of the 
universe; for other limits see Ref. (^j]). 

For the relativistic two-body system of the action-at-a-distance electrodynamics, general 
solutions are not known and the only known special solution is the circular orbit for the 
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attractive two-body problem, first found in Ref. [I| and later rediscovered in Ref. jl(| ( see 
also Q and Q. Our problem na, a, ready been studied: the symmetric motion of two 
electrons along a straight line [—x%(f) = Xi(t) = x(t)], which has the following equation in 
the action-at-a-distance electrodynamics 

d, v f_( \-v{t-r)lc \ e 2 f l+v(t + q)/c \ 

m dt K s J 1 _ {v / c) 2> 2r2 \l + v(t-r)/c) + 2? \l - v(t + q)/c) ' 1 ] 

where v(t) = dx/dt is the velocity of the first electron, of mass m and charge e, and r and q 
are the time-dependent delay and advance, respectively, which are implicitly defined by the 
light-cone conditions 

cr(t) = x(t)+x(t-r), (2) 
cq(t) = x(t)+x(t + q), 

where c is the speed of light. In general, a neutral-delay equation such as Eqs. (JTJ) and (J2J 
requires an initial function as the initial condition, but for the special case of equations (JI} 
and (J2J it was proved in 1979 that for sufficiently low energies the Newtonian initial condition 
[ x(0) = x and v (0) = v a ] determines a unique s ymm etric solution that is globally defined 
(i.e., that does not runaway at some point) This surprising uniqueness theorem 

reducing the initial condition from an arbitrary function to two simple real numbers (initial 
position and velocity) already suggests that the physical phase space could be isomorphic 
to a two-degree-of-freedom Hamiltonian vector field, at least for low velocities (which is 

n 

what we find here). The first numerical method to solve Eqs. (0) and (EI) was given in 15] 
and converged to solutions up to v/c = 0.94. Later another method converged up to 
v/c = 0.99. 

In the following we present a method to find the nonrunaway solution of Eqs. and 
(J2J with a two-degree-of-freedom Hamiltonian system. Our method is based on the physics; 
it starts from the Fokker action and transforms the neutral-advance-delay equation into two 
separate Hamiltonian ordinary differential equations for the same trajectory in two different 
foliations. The Hamiltonians are defined in implicit form and can be solved explicitly in 
terms of the arbitrary ghost functions by using the Hamilton- Jacobi theory. The condition 
that the two solutions describe the same trajectory poses a functional problem involving one 
of the ghost functions, with known asymptotic form, which must be solved self-consistently 
and is the basis of our numerical calculation of the Hamiltonian orbits. We outline the 



3 



method with an emphasis on the physics described by this complex conservative dynamical 
system. The paper is organized as follows: in Sec. II we describe the bi-Lagrangian method 
and solve explicitly for the motion resulting from the Hamiltonians. In Sec. Ill we discuss 
some consequences of symmetry on the explicit solution of Sec. II to reduce the number of 
arbitrary functions. In Sec. IV we determine an equation to match the dynamics of one of 
the particles in the two foliations, which turns out to involve only one of the ghost functions. 
In Sec. V we use steepest-descent quenching to find the self-consistent Hamiltonian orbits. 
In appendix 1 we prove the twice-monotonic property for the arbitrary-mass case. Appendix 
2 discusses an alternative covariant derivation of the equal-mass case, and we also show here 
that the action of a Lorentz transformation on the Hamiltonian is represented by a canonical 
transformation with rescaling of the evolution parameter. In Sec. VI we give the conclusions 
and discussion. 



II. OUTLINE OF THE METHOD 



Here we consider the isolated two-body system with repulsive interaction, away from 
the other charges of the universe, a conservative time-reversible dynamical system in the 
action-at-a-distance electrodynamics. The equations of motion for two bodies in the action- 



at-a-distance electrodynamics of Whee 
Eqs. ((H) and (J2J), are derived formally 
action functional 



er and Feynman |2j, henceforth called 1D-WF2B, our 
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by extremizing the Schwarzschild-Tetrode-Fokker 



S F = -J mi d 3l -J m2 ds 2 -e> J fs^-x^-i^ds,, (3) 

where Xi(si) represents the four-position of particle i = 1,2 parametrized by its arc-length 
Si , double bars stand for the four-vector modulus \ \xi — x 2 \\ 2 = (x\ — x 2 ) ■ {x\ — x 2 ), and 
the dot indicates the Minkowski scalar product of four-vectors with the metric tensor 
(goo = 1, gu = Qw. = <?33 = —1). The particles have masses mi, m 2 , common charge e, and in 
our units, c = 1 The formal conserved energy associated with the Poincare invariance 
of the Fokker Lagrangian (J3J) is discussed in Ref . 0, ^| , a nonlocal expression involving 
an integral over a portion of the trajectory, which is not useful to the present work, even 
though we start from the same Lagrangian (J3j). 
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The starting point of our method is a transformation to new variables 



£x = tx - xt, Ci = h + x i, 

£ 2 = t 2 -X 2 , C2 = t 2 +X 2 . 



(4) 



As first noticed in Ref. 
separate local actions 



171 ]. this transformation splits the action integral (J3j) into two 



(5) 



with 



S a = - J mMidCi) 1 ' 2 ~ J rn 2 (d£ 2 dC 2 ) l/2 (6) 



and 



S b = - J rm^dCi) 1 ' 2 - J m 2 (d^ 2 dC 2 ) l/2 (7) 



2 /■ A 5(6 - 6; 



~ e y y | Cl _ C2 | Ki^ + ^Ci). 

It should be noticed that the double integral of Eq. is reduced, after integration of the S 
function, to a single integral over the parameter of particle 1, with particle 2 contributing only 
at the advanced and retarded positions, this being precisely the reason for the nonlocality 
of the theory, as illustrated in Eq. ((TJ. The usefulness of parametrization (0J) is that it 
naturally breaks the double integral of Q into two integrals, each involving a different 5 
function, and integration over each 5 function couples particle 1 with particle 2 at either 
the advanced position [the double integral included in S a of Eq. ©] or at the retarded 
position [the double integral included in St of Eq. (JJJ)]. For example, in action S a of 
Eq. (jnj), the nonzero contribution of the 5 function occurs where the parameters £i and £ 2 
take equal values £i = ( 2 =(, and this ( is the natural independent parameter of the local 
action S a (this parametrization is often named front form of dynamics in the literature and 
we henceforth call it type a foliation). For action Sb of Eq. (JZJ), integration over the 5 
function produces a contribution to the integral only where the two parameters £1 and £2 
are equal, and £1 = £2 = £ is the natural time parameter of action Sb (henceforth called 
type b foliation). To gain some insights into the two types of foliation, we notice that with 



type a, the particles are automatically in the light-cone condition (x\ —X2) 2 — (£1 — 1 2 ) 2 = 0, 
particle 2 always being ahead of particle 1 in time after the choice x\ — x 2 > 0, with the 
light-cone distance being 

r 8 = ~(6-6)- (8) 

With type b parametrization the particles are also in the light-cone condition, with particle 
2 behind in time and the light-cone distance being 

rfc = |(Ci - Ci). (9) 

The first heuristic guide for this work, as first noticed in Ref. is the simplicity of the 
Euler-Lagrange problem for partial action © : after expressing action @ in terms of the 
time-like parameter (, it is easily verified that the associated Euler-Lagrange equation is a 
simple ordinary differential equation (not a delay equation anymore!). The Euler-Lagrange 
problem for action (JJJ) is analogous, with £ replaced by £. To avoid confusion, we henceforth 
define that a Lagrangian has a local form when the associated Euler-Lagrange problem is 
defined by an ordinary differential equation. In searching for a local Lagrangian problem, we 
could try to extremize each of the partial action functionals of Eqs. (jOJ and (jZJ) and obtain 
a trajectory by solving the Euler-Lagrange equation for either SS a = or SSf, = 0. Each 
separate minimization, in general, yiel ds a different trajectory, which is the time-asymmetric 
problem studied in several works 19]. The main idea of our method is that if these two 
trajectories turn out to be equal, this common trajectory also extremizes the original action 
integral ©, as 6Sf = ^SS a + ^5Sb = + = 0. Simply formulated as above, the problem 
turns out to be impossible; and it is possible to prove that the two separate solutions can 
never describe the same orbit. To overcome this difficulty we need to postulate a more 
general bi-Lagrangian problem by simultaneously solving 

6S a = 6G, (10) 

and 

SS b = -5G, (11) 

with G being a so far undetermined Lagrangian. A trajectory that satisfies Eqs. (fTUjl and 
(fTTj) will also extremize the Fokker action (jHJ), a simple consequence of Eqs. (jSJ, (fTU|) and 
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(HIJ: 

SS F = l -5S a + l -5S b = ^6G - l -5G = 0. (12) 

Our first task is to find a sufficiently general Lagrangian G such that Eqs. (fTU|) and (fTTj) 
yield the same trajectory. Once we are trying to avoid delay equations, it is desirable that 
the Euler-Lagrange Eqs. for (fTT)|) and (fTTf be ordinary differential equations, which is the 
heuristic guide for choosing the functional G. A functional G that leaves the two separate 
problems (jl()|) and (fTT|) in the local form is henceforth called a bilocal ghost Lagrangian 
G. Here we consider symmetric and time-reversible solutions of Eq. only, but for the 
variational calculus that follows, it is necessary to study such an orbit immersed in a family 
of orbits, defined as follows: A time-reversible orbit naturally defines a preferred frame; the 
Lorentz frame where the orbit is time-reversible, and we henceforth call it the center of mass 
frame (CMF). We consider in the CMF the family of all orbits such that the trajectories of 
electrons 1 and 2 are both time-reversible but not necessarily equal (non-symmetric orbits) 
[xi(— t) = Xi(t)] and [x2(—t) = ]> an d with the physical property that both the 

advanced and retarded distances decrease monotonically to a point of minimum and then 
start increasing monotonically again, as illustrated in Fig. 1. We henceforth call this family 
of orbits the CMF family. The fact that the solution of Eqs. (pQ) and (j2J has this piecewise 
monotonic property is a consequence of the velocity being a monotonic function of time, 
which was proved in Ref. [13] for sufficiently low velocity orbits (in Appendix 1 we prove 
this assertion for the arbitrary-mass case). We henceforth refer to a CMF orbit as a twice- 
monotonic orbit. Since the solution we are looking for is symmetric and time-reversible, 
it obviously belongs to the CMF family, and since this solution extremizes Eq. (JHJ) in the 
family of all orbits, it obviously does so restricted to the CMF family. A symmetric and 
time-reversible orbit seen in a Lorentz frame other than the CMF has the property that the 
future of electron 1 is the past of electron 2 and vice-versa. In the following we restrict the 
analysis of the different-mass case to the CMF. For a general covariant derivation of the 
equal-mass case see Appendix 2. 

In the following we prove four integral identities for the orbits of the CMF family, which 
are later used to construct the bilocal ghost Lagrangian. The action of the time-reversal 
operation on orbits of the CMF family can be shown to be the following map: £1,2 — > —£1,2 , 
£i j2 — > — Ci,2j r a — ► r b and it is worth noticing that time reversal maps type a parametrization 
onto type b and vice-versa. In this work we ignore questions of convergence and define all 
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integrals formally from — oo to oo, an ambiguity inherited from the Wheeler-Feynman theory 
and discussed in Ref. jl7|. The simplest type of integral identity we shall use, valid for an 
arbitrary function <p(x) of the real variable, is 

/ <KCR = / 0(O#- (is) 

J a Jb 

The lower index of the integral denotes the parametrization type, and the above identity is 
trivial, as with either type a or type b the parameter runs from — oo to oo (( for type a and 
£ for type b). It is also interesting to look at Eq. (JT5j) as a consequence of the coordinate 
transformation induced by the time-reversal symmetry of the CMF family ( ( — > — £). In the 
same way, we can prove in the CMF the following integral identity, involving an arbitrary 
function V(x) of the real variable 

The combination (d^id( 2 + d&^Ci) * s the time reversible Lorentz-invariant area element 
that appeared naturally in Eqs. © and (f7j). Last, the same time-reversal action (£ 1)2 — ■> 
— £12, £1,2 - ► — 0,2) on the CMF family produces the following identities for arbitrary func- 
tions «(C) and /3(() of the real variable; 

/ a d , t (OKi<iCi) 1/2 = / a«M(0(#i<*Ci) 1/2 , ( 15 ) 

Mcm 2 dc 2 ) 1/2 = [ Pd, t ^M 2 dc 2 ) 1/2 . 



The above identities suggest that we use a ghost Lagrangian G of type 

Notice that the dot over £i >2 in equation (fTH|) indicates the derivative respect to ( ( the 
time-parameter of case a). This G is in the local form when added to S a , where ( plays the 
role of the time parameter and the coordinates are £1, £ 2 . When this same G is subtracted 
from action S&, the integral identities allow us to express G as 

g = j h m) + \v{o{Ci + 6) + «(£)v^+ /9(o VS]^. ( 17 ) 

which is also in the local form for action Sb, with £ being the time parameter and the 
coordinates being d and ( 2 . Notice that the dot over £ 12 in equation (|T7|) indicates deriva- 
tive respect to £ (the time-parameter of case b). One could in principle add more general 
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parametrization-invariant terms to G ; for example, terms involving the integration ele- 
ment (^ic/Ci) 1/4 (^ 2 ^C2) 1/4 

or any highly composite term, and the inversion to the Hamil- 
tonian formalism would involve several branches. Lagrangian is the most general ghost 
Lagrangian whose associated Hamiltonian involves quadratic rational functions of the mo- 
menta, and should suffice if the orbit has only two monotonic branches, corresponding to 
the two elements of the Galois group of a quadratic equation. The need for only four func- 
tions becomes also clearer later on, when we find that there are four determining equations 
involving these four arbitrary functions. We notice also that <fi is defined up to a constant 
in Eqs. (jlfij) and Q17j) . which is also true of V, as adding a constant to V simply adds a 
total time derivative to G (a gauge transformation). There is also no gain in generality if 
one defines a general linear term like (ViCi + V2 C2) m Eq. (|17|) . as this is also a trivial 
transformation of the case we used. 

In the following we guide the reader to a division of the phase space into two disjoint 
regions, as our constructive method defines one Hamiltonian for each separate region as an 
implicit function of phase space: The condition r = divides the phase space of a twice- 
monotonic orbit in two separate regions according to whether r > or r < (in Appendix 
2 we show that this splitting is actually a covariant splitting for the equal- mass case). The 
change from ( to r a is one to one in each of the two regions of a twice-monotonic orbit 
C € [— 00, — |£c|] and ( G [— \( c \, 00], as can be seen from Figure 1, and this naturally splits 
all integrals into two. For example, the left integral of identity |T3|) splits as follows: 

/ <KCR = r' Cc '<Mr a R+ I" <Pt(r a )d(. (18) 

Ja J— 00 J — |f | 

The above integral identity involves two arbitrary functions, and the subscript t (as in turn) 
indicates that function <p t is defined in the region of phase space where the trajectory of 
particle 1 includes a turning point (see Fig. 1), while subscript d (as in direct) indicates 
that <fid is defined in the region of phase space where the trajectory of particle 1 is without 
a turning point (see Fig. 1). The same integral can be expressed for type b applying the 
time-reversal change of variable £ — > — £ to the right side of Eq. (|18jl. which maps r a to r& 
and maps the critical point ( = — |Cc| of the ( parametrization to the critical point £ = \( c \ 
of the £ parametrization, 
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/r\Qc\ roc 
<K£R=/ Mn)<%+ Mn)<%. (19) 

These two portions are indicated in Fig. 1 for both case a and case b. In the following we 
split all integrals of the ghost Lagrangians (|16|) and (JT7j) in two, which we indicate with 
subscripts t and d in the same way as of Eqs. (fT£J) and fTT)j) . The usefulness of the above 
splitting of the phase space is that one can express all functions in the ghost Lagrangian as 
functions of the light-cone distance (JEJ) in each region, such that the Lagrangian becomes 
independent of the time-like parameter and allows the existence of a conserved energy E. 

Our next task is to solve each separate problem for each separate Lagrangian (and for each 
region of the phase space). After inclusion of the ghost Lagrangian, the problem SS a = SG 
in the CMF family implies the Euler-Lagrange equations for L a = S a — G, 



L\ A = -J [Ml^ + M^^+ij^^ + iy w (r.))(6 + 6) + Mr a )]dC, (20) 

where Ml' a = m\ + a tt d{r a ) and = m 2 + A,d( r a) and the Lagrangian can be uniquely 
inverted in each branch to produce a Hamiltonian, because of the monotonic property. The 
problem SSb = —SG is described by Lb = Sb + G, 



L t4 = _J [Mgyfe + Mgyfe + (t^tt - ^ M (r 6 ))(Ci + Ca) - MnM, (21) 

with M l { h d = m\ — a ti d{rb) and M^f = m 2 — (3t,d{ r b) ■ We have introduced eight arbitrary 
ghost functions: 4>t,d , Vt,d > «t,d) and j3 t ,d , four for each separate region of phase space, and 
we notice that these ghost functions enter with a plus sign in case a and with a minus sign 
in case b. The Hamiltonian in each case is given by 

H a = ~Trr7 ~ 2 \ + , ~T#g e2 J - ( 22 ) 

and 

#6 = — { ; — - — 2 — + \ — - — 2 — } + Hn)- 23 

We have omitted the subscripts but it should be kept in mind that each of the above 
Hamiltonians is defined separately in each region of the phase space, a separation that will 
be useful when we come to the symmetry considerations. Notice that the Hamiltonian H a 
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depends only on r a = — — £2), which implies that P a = p\ + p 2 is a constant of motion. 
For type b parametrization, Hamiltonian depends only on = — ( 2 ), implying the 
constant P& = p\ + p 2 . The constant P a = p\ + p 2 suggests a canonical change of variables 
for Hamiltonian H a , defined by 

X = + P = Pi+P2 (24) 

x = -(6 -6), p = pi-P2- 

For type 6 we use the analogous transformation with £ replaced by £ in the above formulas. 
One can use Eq. (J2l|) to express pi and p 2 of Eq. (}22|) in terms of the constant P = P a and 
the relative momentum p, and substitution into the condition H a = E a yields a quadratic 
equation for p, with solutions 



» = ± + i^j + ,/(r) + ^ + 1 (trt J • (25) 

where Q a = j(^i a + -^fa) anc ^ ^ a = \(^2a ~ ^ia) an d r a = M- The separation for case b 
is analogous. 

So far we have shown that any common solution of Hamiltonians f!22l) and (|2Hjl is also 
a solution of the original advance-delay problem of Eqs. (0) and (0), for arbitrarily given 
potentials (p t ,d, V tjC i , ott,d, and It turns out that, even if we guessed the four potentials 
correctly, Hamiltonians (j22j) and ()23|) would have only a single trajectory in common for each 
given set of potentials (this becomes clear in the numerical work of Sec.V). This obstacle 
can be overcome with the Hamiltonian formalism if we generalize the potentials of Eqs. fl22(l 
and (j23|) to implicit functions of the energy E a = H a in case a and of = Hf, in case b. For 
example, the potential is generalized to = <f)(r a , E a ) in case a and to = 0(rb, Ef,) in 
case b (an analogous generalization goes for V, a, and (3). This generalized Hamiltonian is 
still a function of phase space, because E itself is a function of phase space, even though it 
is now only implicitly defined by Eqs. (J22)) and (J23j) . In this generalization, for each given 
orbit, of energy E a , we still define the ghost Lagrangians with (ffllj) and (J2"T|) using fixed form 
potentials: = 4>(r,E ); V = V(r,E ); a = a(r,E ); and (3 = (3(r,E ). By construction, 
these generalized ghost Lagrangians have only a single orbit in common with the generalized 
Hamiltonians, but it is essential that such provisional Lagrangians exist, such that we can 
prove that the Hamiltonian equations associated to ()22)1 and (I23jl reduce to Eq. (pj, which 
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is accomplished by using Eq. (|12|1 with fixed form potentials. After that we can dispose 
of the Lagrangians. On the Hamiltonian side, if we are changing the potentials with E the 
Hamiltonian equations of motion derived from (|22|) and (|23|) pick an extra term proportional 
to the derivative of the Hamiltonian with respect to E (due to the implicit dependence). 
We must therefore supplement a condition that this derivative vanishes along the orbit in 
each case, which in case a reads 

dH a ( Pl P,r,E a ) 



and, in case b 



dE a 
dH b (p } P } r,E b ) 



0, (26) 



0. (27) 



dE b 

In Eq. (I2fjjl . the derivative dH a /dE a = should hold on the energy shell H a = E a and in 
Eq. (j27j) the derivative dH b /dE b = should hold on the energy shell H b = E b . Elimination 
of the relative momentum p a from H a = E yields Eq. 1)2*5)1 and substitution of ([2*5)1 into 
(|26|) yields a partial differential equation (PDE) involving the four potentials. An analogous 
PDE results for ([2*7)1 in case b such that Eqs. ([2*"""l and ([2*7)1 define two partial differential 
equations involving the four arbitrary potentials in each region (variables of the partial 
differential equations are r, P, and E). Rigorously, the generalization to implicitly defined 
Hamiltonians proceeds only if the time-reversal operation also maps E a into E b . For that 
we notice that <fi enters with a plus sign in H a [equation (|22|1] and with a minus sign in 
H b [equation 1)23)1]. and the required symmetry can be accomplished by adding an energy 
dependent constant to 0. We conclude this paragraph stressing that the generalized ghost 
Lagrangians were only a provisional artifact en route to the eventual derivation of the 
Hamiltonians ([2~2*)l and (J2Hj) from a variational argument with use of symmetry. It should 
be clear that after we generalize Eqs. ([2~2*)l and ("2HJ) to implicit dependence and postulate 
Eqs. (J26)) and ([2*7)1 we can no longer go back to the simple provisional ghost Lagrangians, 
and our constructive approach is essentially left with an implicitly-defined bi-Hamiltonian 
system. 

In the following we show that even with Hamiltonian ()22jl defined in the implicit form we 
can write out the motion explicitly: This explicit solution is accomplished in the manner of 
Hamilton- Jacobi, by use of a canonical transformation with a generating function S given 
by 

S = PX + W(x, P, E) - E(, (28) 
12 



where the function W(x, P, E) is defined by integration from the condition p = dW/dx 
, with p given by Eq. (|25j). This canonical transformation is defined such that the new 
momentum associated with the old variable X is the same old constant P = dS/ dX and the 
other new momentum is the energy E (with this last definition we exploit the fact that E 
is already one argument of the potentials). We choose S in the manner of Hamilton- Jacobi 
such that the new Hamiltonian vanishes: K = H + ^| = 0. As the Hamiltonian is zero, the 
new coordinates are defined simply by two constants X and Cq 



X = dS/dP = X + dW/dP (29) 
C = -dS/dE = C - dW/dE 

The above equations for type a define ( and X as functions of the variable r a = \x\, and 
provide the complete solution of the Hamiltonian motion. For further use, it is interesting 
to take the differentials of Eq. (|2T?j) relative to x, 

dX = -(d 2 W/dxdP)dx = -(dpfdP)dx, (30) 
d( = (dW/dxdE)dx = (dp/dE)dx, 

where we have used p = dW/dx (definition of the Hamilton- Jacobi transformation) and 
exchanged the partial derivatives. The explicit form of the differential for the trajectory is 
obtained using Eq. (j^J) to relate particle coordinates to X and ( and using (|3Ujl to relate 
dX and d( to dx. For type a parametrization the explicit solution is 



dtla — 


1 

-(d( a + dX a + dx a ) 


= i,dp a d Pa 

2 [ dP dE > ai 


(31) 


dt2a = 


^{d( a + dX a - dx a ) 


1 dp a dp 

2 [ dP dE + } a ' 




dx la = 


^{d( a - dX a - dx a ) 


_1 ( d Pa dp 
2 [ dP dE } a) 




dx 2a = 


^(d( a - dX a + dx a ) 


l.dp a dp a 

= 2 { dP + dE +1)dr - 





where we have also used dx a = —dr a . Analogously for type b (dxb = dr b ) we obtain the 
explicit solution 

dt lb = ~(d£& + dX b + dx b ) = ~(|| - ^ + l)dr b , (32) 
dt 2 b = 7^(dt,b + dX b - dx b ) = 7^(^| - ~ l)dr b , 
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dXlb = \(~ d & + dX b + dx b) = ^(1 - |P - ^) rfrfe ' 

We recall that Eqs. (|31|) and (J32j) give the explicit solution in terms of p(r, E, P) as given 
byEq. 



III. SYMMETRY CONDITIONS FOR THE EQUAL-MASS CASE 

In this section we discuss only the equal-mass case, and for that we set mi = m2 = m = 1 
and allow only the charge to be arbitrary, from which the general m and c case can be 
recovered by simply replacing e 2 by (e 2 /mc 2 ). We henceforth set e = 1 as well, which can be 
accomplished by a rescaling of distances accompanied by a rescaling of time to keep c = 1 . 
In the following we derive general symmetry relations involving the eight arbitrary functions 
oct,d , Pt,d , Vt,d , and <p t ,d ■ Formula (|2*5Jl for p a is the solution of a quadratic equation and 
defines two different functions p a by taking the plus and minus signs of the square root. It is 
easy to show that at the branch point the square root vanishes, so that a single branch of the 
square root describes each of the t and d physical regions of phase space as indicated in Fig. 
1, henceforth indicated by p\ and p d a (on p a ^ we use superscripts to indicate branch type, 
to avoid overloaded notation, but with the ghost functions we keep using subscripts). We 
have assumed that the orbit is time- reversible in the CMF and, to be consistent with that, 
time reversal must map each branch of the type a trajectory of particle 1 onto a branch of 
its type b trajectory, with the corresponding velocities transforming like vff b (r) = —vf£ a (r) 



for r G [r , oo], as illustrated in Fig. 2. In an analogous way, for particle 2 we should have 

d,t 
V 2a,b 



' ,A = —v^bair) for r £ [r Q , oo]. These two symmetry conditions, when expressed in terms 



of dp/dP and dp/dE using Eqs. (j31|) and ([32)1 . imply the four conditions 

r\ t,d r\ t,d 

d Pa,b 9Pb,a 



OP dP 

t,d a t,d 



(33) 



i t.d o t.a 

d Pa,b OPb,a m x 

dE dE ' 1 } 

Conditions (J3~3j) and (|34j) represent two conditions for region t and two conditions for region 

d, each involving the corresponding set of four potentials. For example, in region t condition 
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is a simple algebraic equation because the potentials do not depend on P explicitly (a 
possible physical choice on the CMF), while condition (J53j) is a partial differential equation. 
We have completed the determining equations for the potentials, which, for example, in 
region t is composed of (J26|) and (|27|) together with the t sector of Eqs. (|33j) and (|34j) . 
At this point we notice another reason to include only four ghost potentials, as we found 
four determining equations to be satisfied [Eqs. (J25J) . (|2T|). (}3"3"j) and ]. The solution to 
these determining partial differential equations should determine the potentials in the CMF. 
This solution is elaborate and still involves arbitrary initial functions of r, which must be 

n 

determined numerically, which we discuss elsewhere [20]. In this work we calculate the 
Hamiltonian orbits directly with an independent numerical method. 

Finally we notice a symmetry relating region d of case a to region t of case b of the 
equal-mass case: As the direction of time in the CMF is arbitrary and the particles are 
identical, the Lagrangian for the d branch of case a must be equal to the Lagrangian for the 
t branch of case b with particles exchanged and vice-versa, which implies 

oc t = ~Pd, Pt = -014, (35) 
4>t = -<f>d, V t = -V d , 

A consequence of Eq. ([55]) is that E\ = E$ and E% = El as well as P l a = P^ and = P£. 
Because and V are arbitrarily defined up to gauge constants, and as <fi and V enter with 
a plus sign in case a and with a minus sign in case b, we can also choose E t a = E\ = E and 
P l a = P£ = P, such that one can use a common value for all the energies and a common 
value for all momenta, throughout the four combinations of region and case. We henceforth 
indicate energies simply by E and momenta by P. 



IV. EQUATION OF MATCHING FOR THE EQUAL-MASS CASE 



In this section we introduce a simpler description in terms of two simple functions s(r, E) 
and F(r, E) that are immediately accessible numerically. For example the d sector of ()33|) 
and is studied by defining dp d a /dP , dp^/dE, dpf/dP, and dpf/dE in terms of s(r, E) 
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(36) 
(37) 



and F(r, E) as 

Qpd dpf _ cosh [s(r, E)\ 
= aP = ~sinh[s(r,£)]' 
dpj dpj _ F(r,E) 
~d~E ~ ~dE = sinh[s(r,£)]' 

For branch t, the consequences of definitions (J3HJ) and (J37|) and the symmetry relations of 

(j35j) are: (i) that branch t involves the same function F(r, E) of branch d and (ii) that 

branch t involves the function s(r, E) of branch d with a change of sign. The intuitive 

picture is that ±s(r,E) and F(r,E) describe both case a and case b, exchanging branches 

in the same case replaces s(r,E) by —s(r,E), while exchanging case for the same branch 

leaves functions s(r,E) and F(r,E) unchanged. In the following we drop the dependence 

on E of the functions for brevity. 

Now we must impose that the same orbit is a solution of both Eqs. (|22)l and ()23j) . which 

demands that the d portion of the a orbit of particle 1 should coincide with a piece of the 

t branch of particle 1 in case b (see Fig. 1). Notice that this is not the one-to-one branch 

correspondence of the symmetry considerations of Sec. Ill, and we stress the word piece, 

because the branches are changed at different points, as can be seen from Fig. 1. We shall 

henceforth drop the subscript notation, and simply write s a and s b , meaning the plus or the 

minus branch of the function s(r), wherever it applies. We can use Eqs. (|3*Tj) , ([3*2"]l and 

(JUnj), (JBTj) to express the differentials of the particle-1 coordinates with type a foliation in 

terms of r a and s a = ±s(r ) : 

dti a + dxi a = . 1 F(r a )dr a , (38) 
smh(s a ) 



dtia - eta i a = — —- — -dr a 
smh(s a ) 



and with type b, 



dt lb + dxi b = eX ^[ Sb \ dr b , (39) 
smh{Sb) 

dt lb - dx u = . —. . F(r b )dr b . 
smh(Sfe) 

At this point it is convenient to introduce still another function: the velocity function of 
particle 1, which must be the same in the corresponding branches of each case. We define 
it in case a by vi a (s a , r ) = tanh($) and in case b by vi b (s b ,r b ) = tanh($), which yields 

exp(2$) = exp(-s a )F(r a ) = exp(s b )F- 1 (r 6 ). (40) 
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The first condition of matching for the trajectory of particle 1 as described by the two 
foliations, requires that the velocities be the same, Vx a (s a , r a ) = Vx b (s b , r b ), resulting in 

exp(s a ) exp(s 6 ) = F(r a )F(r b ), (41) 

a rearrangement of Eq. (|4Uj) . It is important to stress that differently from the symmetry 
conditions, in this condition r a and r b are not equal, but rather for every pair (s a , r a ) we 
should be able to find a pair (s b , r b ) such that (}4~Tj) is satisfied. Fig. 2 illustrates yet another 
symmetry special for the equal-mass case: while particle 1 has a velocity angle $t(r a ) (event 
v\ a in Fig. 2), particle 2 has a velocity angle of $d(r a ) , the same velocity particle 1 had in 
the past at the first time that r a = r [this symmetry reads vl a (r) = vf a (r)]. Notice that 
there are two points along the orbit where the advanced distance assumes a given value, 
one in the t branch where the velocity angle of particle 1 is $t(r a ) and one in the d branch 
with velocity angle $^(r a ), as illustrated in Fig. 2. With the understanding that these two 
branches must be produced with opposite signs for the function s(r), equation (f4~Uj) implies 
that 

exp($ d ) exp($t) = F(r a ), (42) 

which in turn shows that F(r) is determined by past data only, namely the function $d(r a ). 
Another consequence of Eq. P2"J) is that F(oo) = 1, as $d(oo) = — $ t (oo), the asymptotic 
boundary condition on the CMF. Once the orbit is described by two differentials, there is 
another condition for the orbits to be parallel at all times, which is most easily expressed 
by equating the relativistic proper time of particle 1 in the two foliations: 

(dnf = ^$^F(r a )(dr a y = ^^F(r b )(dr b ?. (43) 
smh {s a ) sinn (s b ) 



From the above we can derive differential equations for the motion of r a and r b 

dr a exp(-s a /2) sinh(s a ) 



dr b exp(— Sft/2) sinh(st) 



(44) 



Notice that we have used opposite signs for the evolution of r a and r b , the only sensible 
choice. Equation (j4^|) describes a decrease of r a and r b at large distances if s a > and 
s b < (ingoing asymptotics) and an increase of r a and r b at large distances when s a < 
and s b > (outgoing asymptotics). While asymptotically s a and s b must have opposite 
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signs, they do not change sign at the same point and, in particular, in the turning region of 
particle 1 they have the same sign, as illustrated in Fig. 1. We can also eliminate s a and Sb 
in favor of $ from Eq. (J44j) . resulting in 

To close the dynamical system of matching we need an equation for the variable $; which 
is provided by the Wheeler- Feynman equation of motion To obtain a local equation, 
we write Eq. ((TJ using a combination of type a and type b foliations in the following way: 
whenever we need the advanced position of particle 2, we write it using type a foliation (as 
particle 2 is naturally in the future light-cone), while the retarded position of particle 2 is 
simply written with type b foliation (where particle 2 is naturally in the past light-cone). 
The usefulness of the variable $ is discovered when Eq. is written in terms of $ and the 
proper time of particle 1, which yields simply (recall that we are using e = 1) 

d$ = 1 exp(2$) exp(-2$) 

dr x 2 i rlF\r a ) + r 2 F 2 (r 6 ) 1 ) 

Equations (}4"5j) and (J4T)|) constitute the complete ordinary differential equation (ODE) to 
describe the matching for the orbit of particle 1. By now we have turned equation ((TJ 
upside down and used all the symmetries, and the resulting Eqs. (|4"5j) and (jH|) are much 
simpler to solve than equation Eq. (JTJ). Rigorously speaking we now have a delay-only 
equation, as F(r) depends only on past data via (|4"2*J) . This should be contrasted with 
Eq. (Q), a neutral-delay-advance equation with infinite lags. To solve Eqs. (|45|) and (|46|) 
one needs to postulate an arbitrary positive function F(r) with a given asymptotic form 
F(oo) = 1, and solve the resulting ODE. For self-consistency the ghost function F(r) must 
be chosen such that the orbit of particle 2 is the same even function as that of particle 1, 
by definition of the CMF. This functional problem is solved numerically in the following 
section. 

Last, as an illustration, we outline how one can express the potentials in the CMF in terms 
of the numerically accessible functions s(r, E) and F(s, E) along the simplest analytical 
solution of Eqs. f[27)j). (|27 |) . and (|34~ |) . As the potentials in the CMF do not depend on 
P, Eq. is algebraic and can be solved simply in both case a and case b. One possible 
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solution is 

Pa , y , e 2 _ [2(1 + g)(l + g) cosh( g ) - (1 + a) 2 - (1 + /3) 2 ] ^ 

a r 4(P + 0) 

and 

p y e 2 _ [2(1 - a)(l - /?) cosh(s) - (1 - af - (1 - /3) 2 ] 

For case a the equations of motion for £i and £2 derived from Eqs. (I22p and (j2*Uj) are 

d£i exp(s) (1 + a) 2 

rfC = ^(r) = 4(p 1 + |V + f)2' 

and 



(49) 



d£ 2 exp(-s) (1 + /3) 1 



~ F(r) " 4(p 2 + ±V+g)2' 



(50) 

2 ' 1 2r ) 

We now take the most physically sensible square root of Eqs. (|%9~j) and (JHU|) and substitute 
into Eq. (J2~2~j) . yielding 

Pf = -0 d + _L[(1 + «,) exp( S /2) - (1 + /3 d ) exp(- S /2)]. 

In the same way, for case b, we obtain after a choice of sign for the square root 

E d b = <p d - -L[(l - a d ) exp( S /2) - (1 - fa) exp(-s/2)]. 
2v r 

As discussed below Eq. (j35|) . we henceforth set Ef = E% = E and P d = P d = P, which 
defines the explicit solution for <p d , V d , a d and [3 d as 

sinh(s/2) ^ 



'F 

V d = VFanh(s/2). 
The resulting equations for a d and [3 d are 

a d exp(s/2) - /3 d exp(-s/2) = 2EVF, (52) 

2(P + — ) 

a d exp(-s/2) - /? d exp(s/2) = -=?— , 

v 

where P = ^{E^-\- Ef) and P = |(P^ + P fe d ) are the sole two physically meaningful constants 
of the problem. At the shortest light-cone distance r Q , which happens at s(r ) = , the 
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determinant of the linear system (|52|) vanishes and poses the following condition involving 
E and P 

(P+-) 

E = - - r ° . (53) 

For r > r Q one has s(r) 7^ and (JH2|) can be solved for ad and (3d yielding 

\EF(r) exp( g /2) + (P + g) exp(- S /2)] 

«d = ^= , (54) 

'F(r) sinh(s) 



[EF(r)exp (-s/2) + (P + ^) exp(s/2)] 
Pd — 



F(r) sinh(s) 

In the CMF P is a function of E, so that the potentials given by Eqs. (joTj) and (p)4"j) depend 
only on E and r. We warn the reader that this simplest type of solution of Eqs. ()26j) . (|27j). 
(j33"| and plj) does not correspond to the low-energy orbits of the repulsive case studied 
here numerically in the following section, and was included as an illustration of how s(r, E) 
and F(r, E) can determine the potentials in the simplest possible way. There are several 
other solutions to Eqs. (j2Hj). (|2T|). (|3*3*j) and and the detailed consideration of all cases 
and energy ranges requires further elaborate analytical work and extensive comparison with 
numerics, and will be published elsewhere. 

V. NUMERICAL INTEGRATION: THE STEEPEST-DESCENT METHOD 



In Sec. IV we saw that describing the same particle 1 in both foliations results in Eqs. 
f!45|) and (|46|). involving the single unknown ghost function F(r). The above discussion 
suggests the following simple self-consistent method to obtain the symmetric solution of 
1D-WF2B in the equal-mass case in the CMF: We start by postulating the functional form 
of F(r), which must go to 1 at large distances, as noted below Eq. (J42)) . For the following 
numerical work we use up to 18 arbitrary coefficients to approximate F(r) by a truncated 
power series, 

n=18 u 

no = 1 - e % (55) 

71=1 ' 

which has the desired asymptotic form. After we assume given values for the k n , the main 
equations (J45|) and (jUp) yield a simple initial value ODE problem. The integration can be 
carried out from the turning point, where T\ = 0, $ = and r a = = r c > r a . (the given 
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functional form of F(r) and r c determine all the subsequent dynamics). Notice that this 
integration automatically produces a time- reversible orbit for particle 1: It can be seen by 
inspection of Eqs. (J45)) and (|46|) that replacing T\ by —T\ and $ by — $ simply exchanges 
r a and r&, such that r a (— tl) = r&(Ti), a consequence of the symmetry imposed. When we 
start particle 1 at the turning point, $ = 0, particle 2 is described in type a parametrization 
at the advanced point (x2 a ,^2a) = ( — f c ,r^), while with type b parametrization particle 2 is 
at the retarded position (x2 a ,^2a) = {—r c , — r c ). As illustrated in Fig. 3, for a generic choice 
of F(r) the future of type b trajectory of particle 2 will not match the type a trajectory of 
the same particle 2, which starts ahead of case b, and the scheme produces two different 
orbits for particle 2, which is an absurd. It is necessary to adjust the function F(r) precisely 
to obtain a single trajectory for particle 2, and it is nice to observe that the asymptotic 
condition F(oo) = 1 guarantees the asymptotic velocity of particle 2 to be the same in both 
foliations, so if we adjust the orbits to overlap in the turning region, they become close 
everywhere. 

Our numerical method produces two trajectories for particle 2 from each set of k n , by 
direct numerical integration of the main equations (|45jl and (|4fi|) accompanied by the driven 
equations for the trajectories: (^fr-) , ( ^frf ) , ) and (^) as determined by Eqs. ([3T]) 
and (J55j). We calculate the trajectories numerically by using a 9/8 embedded Runge-Kutta 
pair. In general, two different trajectories are obtained for particle 2, as illustrated in Fig. 
3 and we calculate numerically the average squared deviation of the two trajectories over a 
grid of positions 

A(k) 



1 N 



\\ Jf 52(ha(x2i) ~ t 2 b(x 2 i)) 2 . (56) 

After that we implement a steepest-descent search in the 18- dimensional parameter space 
governed by the quenching equation dk n /ds = —dA/dk n until it finds a minimum value 
for the squared deviation of Eq. (J56j) (see Ref. 2l| for an analogous numerical quenching 
procedure). 

In Fig. 4 we put the converged trajectory of particle 1 superposed to the reflected 
trajectory of particle 2, for velocities v/c = 0.46, v/c = 0.54 and v/c = 0.71. Notice 
that the trajectories coincide perfectly, indicating that the quenching search satisfied all the 
symmetry relations and thus found a consistent solution. As the solution is self-consistent, 
we can not set the asymptotic velocity directly, and we determine a final low velocity by 
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using an initial condition with r c 3> 1, while a large asymptotic velocity is achieved by 
using r c ~ 1 (one classical electronic radius). We start the numerical work at low velocities, 
by setting a large value of r c , which results in a small asymptotic velocity. After that, we 
decrease the value of r c and give the formerly determined solution as seed to the quenching 
method, which converges much faster as there must be a solution in the neighborhood of 
the seed. 

In Table 1 we list the final velocities v/c as a function of the initial condition r c . From 
the numerically converged trajectories we calculate the minimum distance r and Table 1 
shows r a /r c . We observe that some coefficients k n converge to a value below the numerical 
precision of 10~ 14 such that only the first iV coefficients are significant to the numerical 
precision. This number increases with the asymptotic velocity, as can be seen from Table 1. 

In Fig. 5 we plot the trajectory of v/c = 0.80 , which took twelve hours of numerical 
quenching to converge and still one can observe a slight mismatch of the orbits in the turning 
region, indicating the slowness of the convergence process. For this case, the numerical 19th 
coefficient is still important, indicating that our basis is failing to converge to the solution, 
which is suggestive that something physical is happening above v/c = 0.71, maybe we are 
even loosing the twice-monotonic property (according to Appendix 1 this could happen at 
any point above v/c = 0.33). There might also be another trajectory in the neighborhood, 
which interferes with the convergence, and last, at high energies the form of F(r) as given 
by Eq. ( l55|) becomes too singular at the collision and we need a regularized numerical 
method. Further numerical studies with a regularized numerical method are needed to 
determine if some special bifurcation is happening to the orbit above v/c = 0.71. In this 
high-velocity region the functional problem posed on F(r) might not have a unique solution, 
and for relativistic velocities it is likely that it does not. The symmetric solution was actually 
proved to be unique only up to a small velocity and we managed to go much above the 
low limit set by Driver in Ref. jlflj . In Fig. 6 we plot F(r) versus (r D /r) for the asymptotic 
velocities v/c = 0.46, v/c = 0.54 and v/c = 0.80. Notice that the functional form of F(r) is 
approximately a linear function of (r Q / r) at low velocities, but at larger velocities it becomes 
highly convoluted. 

As a test for the existence of other types of orbits, we integrated Eqs. (j43j) and (j4^j) using 
two completely general (and possibly different) functions F a (r a ) and i^r?,), each defined by 
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an independent truncated power series like 

n=9 t,a,b 

FaAr) = K* ~ £ (57) 

n=l 1 

Notice that the saturation value is not anymore k^ ,b = 1 like in Eq. (|55|) . but rather a 
generic quenchable value in each case. The integration procedure was started with r a ^ rf, 
at $ = for particle 1. We found that the quenching method converged to functions F a (r a ) 
and i*&(r&) that were generally different but were always related by the scaling discussed in 
Appendix 2 [see equation (}65|) and text below it]. According to the discussion of Appendix 
2 this is the case when the orbit is symmetric in another Lorentz frame. This above-defined 
search would be capable of finding any existing orbit with the twice-monotonic property and 
the fact that it always converged to Lorentz-transformed symmetric orbits is indicative that 
there are no other types of low-energy solutions. 

At this point it is interesting to appreciate the big detour taken by our numerical method 
to solve Eq. which should be compared to the most straightforward way to solve a 
neutral-delay differential equation like Eq. ( | T |) . namely: postulating an initial function and 
continuing the solution by use of the differential equation. The straightforward method nec- 
essarily leads to runaways because one is never capable of guessing the unique nonrunaway 
initial function, and even if one does guess the nonrunaway condition, numerical roundoff 
plagues the integration and one still gets runaways after some short time. Our numerical 
method is superior in this respect precisely because it is already placed in the nonrunaway 
manifold and the quenching implemented to solve the functional problem for F(r) is numer- 
ically stable, as it does not involve extrapolation. 



VI. CONCLUSIONS AND DISCUSSION 



In this paper we discussed the solution and Hamiltonian description of the time-symmetric 
two-body problem of the action-at-a-distance electrodynamics with repulsive interaction. 
Our method is closed and does not involve expansions, only the hypothesis that the orbit is 
twice-monotonic was used. We conjecture that our solution is already the general solution 
at low energies, which can be argued in the following way: For a generic solution, possi- 
bly unsymmetrical, one can always find a Lorentz frame where the asymptotic outbound 
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velocities of the two particles are opposite (the outbound CMF). Now if asymptotic data 
determine a unique trajectory for low energies, that trajectory is our symmetric solution in 
this CMF. Of course the symmetric solution is always a possible solution, but the symmetry 
of the physics in this outbound CMF suggests that it is the only solution for low energies, 
as the solution has to correspond to the Coulombian solution, which has this property (this 
would actually be a non-trivial generalization of the work of Driver and Hoag in references 

13 , 1^| )• With the above in mind, the Hamiltonian we derive here is already the general 
order-reduced Hamiltonian for low energies. Some few numerical experiments have suggested 

ihat the conjecture is correct at low energies. Even the high energy solutions found in Ref. 

r| exhibit the property that the future of one particle is the past of the other, and it looks 
like the numerical method in Ref. [l6[ simply picked a Lorentz frame slightly off the CMF, 
but that requires further investigation with a regularized numerical method. 



The idea to remove the field degrees of freedom goes back to Dirac |7j and later Wheeler 
and Feynman planned to quantize WF2B as a means to avoid the divergencies of QED, as 
in the action-at-a-distance theory the infinite number of field degrees of freedom is absent. 
History says that the famous seminar that never came from Wheeler (see Ref. page 
97) was due to difficulties in converting the Fokker Lagrangian (jSJ) to the Hamiltonian form. 
This task is still not fully done and in this work we took a step in that direction for the 
one-dimensional case at low energies. Notice that the implicit dependence of the Hamil- 
tonian operator is actually convenient for an eigenvalue equation, and one could discuss a 
canonical quantization procedure based on either Eq. (j22J) or Eq. pHjl . using the numer- 
ically determined potentials. Of course Wheeler and Feynman were mainly interested in 
the attractive case, of greater relevance for atomic physics and specially for the Lamb shift 
calculation. The attractive problem is being published elsewhere j^j]. In this same chapter 
5, page 97 of reference zj, Feynman says that ' I didn 't solve it either — a quantum theory 
of half- advanced half-retarded potentials — and I worked on it for years... '. This is still an 
outstanding problem today and the difficulties in casting relativistic Lagrangian interac- 
tions into the Hamiltonian form are well explained in references The only studies 
we know of dealing with the time-symmetric problem involve power expansions. We are 
aware of another attempt at a Hamiltonian description of 1D-WF2B that ends up with an 
infinite-dimensional Hamiltonian 2^|, such that further order reduction is needed to select 
nonrunaway orbits. 
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Our description might seem to violate the no-interaction theorem , but there are 
two places where we avoid it: (i) the no-interaction theorem is an obstacle to covariant 
Hamiltonian description of two interacting particles only in fully 3-dimensional motion. We 
are restricted to 1-dimensional motion; (ii) The evolution parameter in our Hamiltonians is 
not time but rather £ for case a or £ for case b and, therefore, the no-interaction theorem 
does not apply. In principle, because time is not the evolution parameter, even in three 
dimensions the no-interaction theorem would not be an obstacle to an analogous procedure, 
and that is an open problem. 

As regards the applied mathematics literature of delay, the theory of delayed functional 



equations 



30[| is a difficult and poorly investigated subject but it turns out that there 



are already a few results worth noticing: In a paper of 1974 by Kaplan and Yorke 
it was noticed that for some special types of delay equations, solutions can be found by 
searching the periodic orbits of an associated ordinary differential equation. This was further 
generalized in 1999 [^J and it was shown that for a large class of delay equations the 
associated ODE turned out to be a Hamiltonian ODE, quite a curious result 2j| brought 
up by applied mathematicians with no relation to either quantum mechanics or Wheeler- 
Feynman electrodynamics. Another set of studies of applied mathematics focuses on the 
similarities of delay equations to either ODE's or extended systems j^j: if the delay is 
small and bounded, the behavior should be reminiscent of that of ODE's, as determined 
by the dimension of the attractor in several systems with small delays, a generic class that 
contains the bound states of the attractive 1D-WF2B, apart from the fact that our system 
is conservative. In the limit where the delays are very large, delay equations are found to 
behave like extended systems, with large dimensional attractors, which is the generic class 
of the repulsive case of 1D-WF2B, where the lags are unbounded and also of the unbound 
states of the attractive case. 

As we mentioned in the Introduction, the action-at-a-distance electrodynamics is capable 
of describing the whole of classical electrodynamics as a limiting case, and even better, 
a limiting case without the complications of mass renormalization, as demonstrated by 
Wheeler and Feynman |3[ . This was actually what led Wheeler and Feynman to the action- 
at-a-distance electrodynamics in the first place, but in doing that they formulated a very 
complex conservative physical theory (the conserved energy associated with the Fokker 
Lagrangian is discussed, for example, in Ref. Q]). It is important to stress that the converse 
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of the above statement is not true at all: the complex conservative dynamics of the action-at- 
a-distance theory is not reducible to a limiting case of Maxwell's electrodynamics (which is 
always a dissipative theory because of the radiation). Relativistic action-at-a-distance shares 
the conservative character with Newtonian gravitation, and in the presence of a universe of 
particles, an atom described by the action-at-a-distance theory has the possibility to behave 
in a way analogous to the solar system in the Newtonian sky: distant solar masses being 
just small perturbations, as opposed to the description a la Maxwell, where it has to radiate. 
It appears to us that the analysis of the complex conservative dynamics of WF2B is bound 
to reveal interesting new physical insights. As we have seen with this special case study, the 
physical nonrunaway condition performs the magical reduction from the infinite-dimensional 
dynamical system posed by the delay equation to a finite-dimensional one, and the large 
body of existing understanding on qualitative behavior of finite-dimensional vector fields 
should be applicable Some results already published for systems of atomic physics 
within the Darwin approximation, a low-velocity Hamiltonian approximation to action-at- 
a-distance, have already revealed interesting news 8ll. Is^ . Existing numerical methods for 
the relativistic case are still short reaching (l2j and can not reproduce the massive numerical 
search performed with the Darwin approximation in 



VII. APPENDIX 1: PROOF OF THE TWICE-MONOTONIC PROPERTY 



In this section we show that in the Coulombian limit of a low-energy orbit, the solution 
of Eqs. and (0) has only two branches, one defined by f > and another defined by 
r < 0. Because of time-reversal symmetry, the theorem is the same for either r a [indicated 
by q in Eqs. (JTJ) and (J2J)] or by r b [indicated by r in Eqs. (JTJ) and (JSJ)]. It suffices to prove 
that there is only one point where a vanishes, with q defined in Eq. A special version 

of this proof was given in Ref. along symmetric orbits of the equal-mass case. The 
proof is trivial and can be done for a generic orbit of the arbitrary-mass repulsive two-body 
system in the CMF: We start from the definition of the light-cone condition for a generic 
CMF orbit, 

q = Xl {t) -x 2 {t + q), (58) 
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where x\(t) represents the position of particle 1, assumed on the right, and x 2 (£) represents 
the position of particle 2, assumed on the left, and we have set c = 1 . Notice that Eq. (J2J) 
is a special case of Eq. ()58j) for symmetric orbits of the equal mass case [x 2 (t) = —Xi(t) ]. 
The key observation is that because the interaction is always repulsive, the velocity Vi(t) of 
particle 1 is a monotonically increasing function of time (particle 1 is repulsed to the right), 
while the velocity v 2 (t) of particle 2 is a monotonically decreasing function of time (particle 
2 is repelled to the left). If we take the derivative of Eq. (|58|l respect to t and isolate q we 
obtain 

• Mt) ~ V2(t + q) 

q = l + v 2 (t + q) • (59) 
For low-energy we have the bounds \v 2 (t + q)\ < v 2 (oo) <C 1 and |t>i(t)| < 1*1(00) <C 1 , 

and therefore the denominator of Eq. (JoTJj) is always positive. In the CMF the value of q 

changes sign from the inbound asymptotic region to the outbound asymptotic region, with 

values onto the interval 

_ v 1 (oo) + v 2 ^) v 1 (oo)+v 2 (oo) 

v l + f 2 (oo) ; v l-u 2 (oo) ' v ; 

To complete the proof we need only to notice that the sum of two monotonically increasing 

functions [t>i(t) and —v 2 (t + q) ] is also monotonically increasing and, therefore, can only 

vanish once. It should be noticed that —v 2 (t+q) is not necessarily a monotonically increasing 

function of t for high-velocity orbits, as 

*^ = (1+i) w (61) 

which is the product of (1 + q) times a negative number [recall that v 2 (t) is a monotonically 
decreasing function of its argument and, therefore, —v 2 (t) is a monotonically increasing 
function]. It can be seen that Eq. (J61|) guarantees that v 2 (t + q) is an increasing function of 
t if (1 + q) is positive, which is the case for a low-velocity orbit. A simple estimate for the 
first velocity where the twice-monotonic property can fail in the equal-mass case is given by 
setting \q\ = 1 in Eq. (joTIj) . which predicts v(oo) = 1/3 . 

VIII. APPENDIX 2: CO VARIANT DEFINITION OF THE EQUAL-MASS CASE 



In this appendix we exhibit a covariant derivation of the above Hamiltonization procedure, 
which we develop only for the equal- mass case and with the hypothesis that the orbit is twice- 
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monotonic (proved in Ref. |l3( for low energies and in Appendix 1 for the arbitrary-mass 
case in the CMF). The definition of a covariant family starts from the observation that the 
Lorentz transformation of a symmetric orbit has the property that the future of particle 1 is 
the past of particle 2 and vice versa. We define the relativistic symmetric family of orbits, 
henceforth called RSF, as the family of orbits with the property that if + £2(^2) — 0, 

then it follows that tx+t% = 0. It is easy to verify that RSF is a Lorentz invariant family of 
orbits and also that any Lorentz transformation of a symmetric orbit belongs to RSF (but, in 
principle, these should not exhaust the RSF: there could be other types of orbits). A generic 
orbit of the RSF is represented in Fig. 7, where we illustrate the time-reversal and exchange 
symmetry. The above definition implies that the future of particle 1 is the past of particle 2 
in the RSF. Inside the RSF, by use of the time- reversal operation £i — > — £ 2 , C2 — * — £1, we 
can prove the same identities (|I3|L ([Bj) and an equivalent form of Eq. (fT3j) . relating particle 
1 to particle 2, (only that in case b the role of a and j3 is exchanged), and these in turn lead 
to the same type of general ghost Lagrangian to describe a twice-monotonic orbit. 

The relativistic condition that the future of one particle is the past of the other implies 
that the solution of WF2B inside the RSF must have the following branch correspondences: 
v ia ( r a) = v 2a( r a) and Ui!?( r &) = v 2b( r b) ( see Fig- These conditions can be seen with the 
help of (}3*T]) and (J3*2|) to be equivalent to the four conditions 



The above conditions imply that the orbit is defined by four different functions ±s (r), 
±Sb,F a (r) and F^ir) [see definition and (J3*7)l ]. To verify that condition (|52Jl is rela- 
tivistically invariant, let us suppose that we tried to describe the orbit from another Lorentz 
frame, with boost parameter to. If the orbit is twice-monotonic, it can be described in 
case a with a Lagrangian L a of the same type of (j2*U|) , and a Hamiltonian of type (|22|h 
and the new coordinates xi a ,i\ a , ii& must be obtained by a simple Lorentz transfor- 
mation with boost parameter w. Imposing this condition on the explicit solution (|31|). and 
noticing that the advanced light-cone distance in the new frame relates to the old one by 




o t,d 



r> d,t 

9Pa,b 



(62) 




dp a 



dE 



(1 - W) dp a 

(1 + w) dE 



(63) 



28 



dp a dp a 

dP ~ dP' 

which, besides showing that Eq. (|62|) is frame independent , also shows that s a (r) as defined 
by Eq. (JSfiJl is a Lorentz scalar: 

s a (r a ) = s a (X a f a ), (64) 
and that F a (r a ) [as defined by Eq. (pITj) ] transforms like 

F a (f a ) = X 2 a F(X a f a ), (65) 



with A a = yj^;- Case b transforms in the same way with A& = 1/A a . This last equation 
allows us to express the Hamiltonian in any frame by use of the CMF form of F(r) and a 
rescaling depending on the boost parameter, which can be determined by asymptotic data. 

Last we show that the action of a Lorentz transformation on Hamiltonian (}2~2*j) is a 
canonical transformation: It is easy to verify with the help of Eq. ( 0] ) that a Lorentz 
transformation simply rescales the coordinates £i and £2 to £1 = j^and £2 — ^£2 with 



A = y YT^' ^ com pl e te the change with a canonical transformation one must scale the 
momenta with the inverse factor, pi = Xp\ and p2 = Xp2- By this canonical transformation 
the transformed Hamiltonian is 

fr ~A r M 2 la {f a ) Mj(f a ) - 

where 0(f a ) = i0(r a ), V{f a ) = XV{r a ) , Mf a (f a ) = Mf a (r a ), and M| (f a ) = M| a (r a ) . 
Notice that Hamiltonian H a picked a multiplicative factor of A and if we also perform a 
change to the natural evolution parameter ( = X( of the new Lorentz frame, it compensates 
exactly for that factor, going back to the form (|22jl , the same form for all Lorentz frames. 
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IX. FIGURE CAPTIONS 



Fig. 1: Particle trajectories in the CMF for ni\ 7^ rri2 in the (x,t) plane, arbitrary units 
with c = 1 and e 2 = mc 2 . Particle 1: case a indicated in branches: la*(solid inner line 
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on right) and la d (hatched inner line on right); case b indicated in branches 16 d (hatched 
outer line on right) and lb 1 (solid outer line on right). For clarity we indicate the type a 
and type b orbits of particle 1 as separate curves, but there is just one orbit for particle 1. 
Trajectory of particle 2 is represented by the solid line on left, branches are not indicated. 
The geometric distance in the (x, t) plane between two points in light-cone condition is \[2 
r, with r being the spacial distance, we dropped the \/2 factor and indicated simply r. 

Fig. 2: Particle trajectories in the CMF for m x = m 2 in the (x, t) plane, arbitrary units 
with c = 1 and e 2 = mc 2 . Indicated is the minimum distance r Q and the velocity of particle 
1 in each branch corresponding to the same distance r = r d = r* > r Q . Event vf a is the 
time-reversed point of event vf b and, because of that, r d = r d = r. A special symmetry of 
the equal- mass case: event v\ a is obtained by particle exchange and time reversal from vf a , 
which implies r* = r d = r. All branches of the trajectory of particle 1 are indicated, and 
omitted for particle 2. Notice that we indicated the spacial distance in light-cone, which is 
\/2 times the geometric distance in the (x, t) plane. 

Fig. 3: Matching the trajectory of particle 1 for a generic F(r) : Particle 2 starts from 
the past in case b and, unless F(r) is chosen precisely, its future does not coincide with the 
trajectory of particle 2 of case a. Arbitrary units with c = 1 and e 2 = mc 2 . 

Fig. 4: Numerically determined trajectories in the CMF for mi = m 2 in the (x, t) plane; 
units with c = 1 and e 2 = mc 2 . Three different symmetric trajectories found by the steepest- 
descent method; the orbit of particle 2 is reflected and superposed onto that of particle 1 
to show the agreement: v/c = 0.46, r c = 4.7; v/c = 0.54, r c = 3.7 ; and v/c = 0.71 and 
r c = 2.97 . 

Fig. 5: Numerically determined trajectories in the CMF for mi = m 2 in the (x, t) plane; 
units with c = 1 and e 2 = mc 2 . The orbit of particle 2 is reflected and superposed to that 
of particle 1, for a highly relativistic case of v/c = 0.8 . Notice the slight mismatch of the 
two orbits, due to failure of convergence of the series. 

Fig. 6: Numerically determined functions F(r), rescaled as a function of (r D /r); both 
quantities are dimensionless. Notice that at low velocities F(r) is well approximated by the 
first two terms of the series, F(r) ~ 1 — ^ , but at larger velocities F(r) becomes highly 
convoluted. 

Fig. 7: A typical trajectory of the RSF in a generic Lorentz frame for the equal mass 
case in the (x, t) plane; units with c = 1 and e 2 = mc 2 . The only symmetry is that the 
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future of particle 1 is the past of particle 2. 

Table 1: Numerically calculated asymptotic velocities (v/c), minimal radii (r /r c ) and 
number (JV) of significant terms of the F(r) series as a function of the critical initial distance 
r c at $ = 0. Units of energy are set by c = 1 and e 2 = mc 2 . 
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